Breakdown of an Electric-Field Driven System: 
a Mapping to a Quantum Walk 
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Quantum transport properties of electron systems driven by strong electric fields are studied by map- 
ping the Landau-Zener transition dynamics to a quantum walk on a semi-infinite one-dimensional 
lattice with a reflecting boundary, where the sites correspond to energy levels and the boundary the 
ground state. Quantum interference induces a distribution localized around the ground state, and 
a derealization transition occurs when the electric field is increased, which describes the dielectric 
breakdown in the original electron system. 
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Dynamics of quantum statistical systems driven by 
finite external fields has attracted much attention as 
a typical class of non-equilibrium phenomenon. One 
problem is the "dissipation" arising in electron systems 
driven out of their ground states by strong electric fields 
as studied by many authors 0, llll IE In 
these references the electric field is expressed (via Fara- 
day's electro-magnetic induction) as a time-dependent 
Aharonov-Bohm(AB) flux, and this induces, for strong 
fields, intcrlcvcl Landau-Zener transitions. There, an 
issue is whether the bunch of energy levels around the 
one having the main amplitude can play a role akin to 
dissipation^. The purpose of the present study is to 
propose a mapping of the system onto a quantum algo- 
rism model for studying the problem. We start by noting 
that the driven quantum system and the quantum algo- 
rism model have in fact similarities: The former deals 
with quantum transitions among macroscopic number of 
energy levels while the latter describes the dynamics of 
many qubits governed by quantum logic gates, for which 
powerful analytical techniques are being developed. 
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where cj , 



creates an electron at i = 1,2, ... L with spin a =|, j, 
the AB-flux measured in units of the flux quantum 
<f>(t) = FLt/h represents the electric field F, and V 
the external potential or electron-electron interaction. 
Disordered mesoscopic systems 0, have been 
treated with this Hamiltonian, which has been extended 
to a strongly correlated electron system by three of the 
present authors 0. The adiabatic energy levels against 
(j> have many level anti-crossings (as schematically 
shown in Fig. 1(a)), which come from the disorder 
in disordered (one-body) systems, or from Umklapp 
processes in correlated systems. The system obeys 
the time-dependent Schrddinger's equation starting 
from the ground state |^o) of H(0), where its formal 
solution is fe-'ti ff C^( s )) ds |* ). When the field is finite, 
non-adiabatic Landau-Zener tunneling processes occur 
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FIG. 1: (a) Idealized energy levels of an electron system on a 
ring plotted against the AB-flux <f> — FLt/h, which increases 
linearly with time after the electric field F is turned on at 
t = 0. Two paths (A, B) for the Landau-Zener transition 
among neighboring levels are shown, (b) A mapping to a 
quantum-walk, where the energy levels are mapped to sites of 
the qubits and the ground state to a reflecting boundary. 



at the level anti-crossings, which take place first from 
the ground state and then among higher excited states. 
Authors of introduced a transfer matrix repre- 
sentation to mimic the evolution fe -4 ^ H(<P( s )) ds ^ ^gj-g 
a set of 2 x 2 unitary matrices represent the transitions 
among neighboring energy levels. Extensive studies 
on this model have found the existence of a dynamical 
localization^, or an Anderson localization in energy 
spaceQ. However, the effect of quantum interference, 
which should be the key to the Anderson localization, 
has been studied only numerically 0,0, 00. 

Now we want to point out that the Landau-Zener tran- 
sitions in multi-level systems presented above can be 
mapped to a quantum walk on a lattice. We shall show 
that the localization (in energy axis) of wave functions 
can be studied in terms of an exact solution for the transi- 
tion amplitudes in the quantum walk. Quantum walk is a 
quantum counterpart of the random walk flTl fTlj ) . The 
mapping we conceive here between the quantum walk and 
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the Landau-Zener dynamics is straightforward: A qubit 
on site n > labeled by L, R corresponds to the states at 
the n-th anticrossing point, which moves down (L; with a 
left-going current) or up (R; right-going) in energy after 
the tunneling event (Fig^). One important point we note 
here is that the mapped quantum walk has a reflecting 
boundary, since we cannot walk below the grou nd state in 
energy. In previous quantum-walk studies(12l ll 'l Il4| an 
absorbing boundary was considered, for which the gen- 
erating functions were obtained [Tij. Here we first ob- 
tain the generating function for the quantum walk with 
a reflecting boundary. We have found the existence, and 
an analytic form for, the amplitude localized around the 
boundary (the ground state in the original problem). The 
solution exhibits an asymmetry between the L and R 
states, which represents a finite total momentum. When 
the electric field exceeds a critical value, a dclocalization 
transition (on energy axis) is observed, which we identify 
here to describe the dielectric breakdown in the original 
electron model. 

Mapping — Let the wave function for the n-th en- 
ergy level at time r (measured in units of At/ 2 = 
h/2FL) be *(n,r) = r), ip R (n, r)), where 

i>L(n,T)(ipR{n,T)) has a left- (right-) momentum. Each 
energy level is subject to a Landau-Zener tunneling (with 
certain probability and phase) to neighboring levels in a 
time period At/2. While transitions among more than 
three levels exist in principle, here we restrict ourselves 
to transitions among neighboring levels, for which the 
quantum tunneling can be described by a set of 2 x 2 
unitary matrices [J, 0, y, Q|- The diagonal elements rep- 
resent Landau-Zener tunneling, while the off-diagonal 
ones the backward scattering process. Here we sim- 
plify the problem by assuming that the matrices are the 

( a ^ ) except for the one at the bound- 



ary, U 



a b 
c d 



As we shall see, an interesting struc- 



ture arises in the overall shape of the distribution in the 
bounded quantum walk even within this simplification. Q 
The wave function evolves deterministically following a 
recursion formula, 

*(n,r + l)=P*(n + l,r) + <9*(n-l,r) (1) 
for excited (n > 2) states, while we put 



*(l,r+ I) = P*(2,t) 
*(0,t+1) = P*(l,r) 



Q*(0,r), 



(2) 



between the ground and the first excited levels. Here 
we have decomposed U into P = (qq), Q = (™) and 
Q = (°°-). Equations ©,© define the mapping to a one- 
dimensional quantum walk on a semi-infinite space with 
a reflecting boundary at n = 0. We take the initial wave 
function \*S? ) to be the ground state, i.e., ^(O, 0) = *(1, 0) 
and *(n, 0) = for n > 1. 

Generating function for the quantum walk — Now we 
state our main result for the quantum walk: We can ob- 



tain the generating function for the wave function in the 
reflecting boundary condition, where the generating func- 
tion is defined as ^ L R (0 — ► n;z) = Y^Lq *pL,n(n, r) z r , 
with the summation taken over all the paths starting 
from and ending at n. The obtained expression is 



* L (0 -> n-z) 



4^(0 ->n;z) = 
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/ d\ + (z^ N 



V a 



l-cA r (0^0;z)' 



A+(z) 



Az 2 + 1 - v /A 2 z 4 + 2A(l-2|a| 2 )z 2 + l 



2dz 



(3) 



Here A = detU = ad -be, and A r (Q -> 0; z) = \d\+(z) - 
Az]z/c is the generating function for the quantum walk 
with an absorbing boundary, governed by Eas. 111211 with 
Q = 0|0|. 

In order to obtain the generating function J3J, we first 
express the transition amplitudes S & (0 — > n;r) as a sum 
of the contributions from paths starting from the 0-th 
level (ground state) ending at the n-th level after r steps. 
For example, E b (0 -> 2; 4) = QPQQ + PQ 2 Q + QQPQ 
consists of three paths, where the quantum nature of 
the walk appears as the matrices being complex, i.e., 
different paths interfere with each other. The transi- 
tion amplitude relates the initial vector ^(0,0) to the 
wave function via ^(n,r) = S fa (0 — > n;r)\&(0, 0). We 
also introduce S a (0 — * n; r) as the transition amplitude 
in the quantum walk with an absorbing boundary; e.g., 
S°(0 -> 2; 4) = QPQQ + PQ 2 Q. 

We obtain the transition amplitudes by employing a 
technique called the PQRS method^]- Since the four 
matrices P,Q,R= (gg) and S = (° a ° b ), with P, Q, R, S 
defined similarly, span an orthonormal basis in M2(C) 
with an inner product (A\B) = Ti(A'B), we can express 
the transition amplitude S b (0 — ► n; r) as a linear combi- 
nation, 



E b (0 



-) =^(0^n;r)Q + 6 f (0^n;r)^ (4) 



with the tilde basis, where P and S terms are ab- 
sent. Thus, eq.(@J defines a set of complex numbers 
6«(0 -> n;r) and b f (0 -> n;r), and we denote S«' f (0 -> 
n; z) = J27 3 =o° 9 ' r (^ ~ * n;r)z T as their generating func- 
tions. The generating functions A^ ,r (0 — > n;z) for the 
absorbing boundary problem is defined similarly from 
S a (0^n;r) = a«(0 -> n; t)Q + d f (0 n; t)R. We note 
that for n = the Q terms vanish, i.e., B^(0 — > 0; z) = 
A q (0 — > 0; z) = 0. For a special case U = U, we drop all 
the tilde marks S 6 (0 n;r), B^ r {0 -> n;z), ^- r (0 -> 
n;r) . . .. 

The PQRS method establishes a connection between 
the generating functions via multiplication rules|l5j for 
the PQRS matrices. From a recursion relation between 
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E b (0 -> 0; r) and E b (0 -»■ n - 1; r - 1), wc obtain 

B«(0^n;z) = 1 + cB f (0 -> 0; z) dzB 9 (0 n . - 1; 

fT(0->ra;z) = [l + aET(0 -> 0;z)l dzB r (0 ^ ?i - 1; z). 

With a repeated application of the no-tilde ver- 
sion of Eq.{HJ we obtain B q (Q — > n;z) = 
[dX + (z)/a] n /d 7 B r {0 -► n;z) = [dA+(z)/a]"[(A + (z) - 
az)/(acz)]. Again using (J5J with a relation between 
£T(0 — * 0; z) and A r (Q — > 0; z), we arrive at eqs.© after 
a bit of algebra. 

Edge states in energy space — One asset of the gen- 
crating function is that we can show that there exist 
"edge states" in the bounded quantum walk. To actu- 
ally show this let us take the Landau-Zener tunneling 
probability, p = exp(— irF/F), where the Zener thresh- 
old electric field F <x A 2 depends on the size of the gap 
A [H[l7l. We plug this in the transfer matrix, 



U 



03 



-pe~ 



y/pe-*! 3 1 ' U 



e % ~< 
-e~ iX i 



( 6 ) 

where the diagonal elements of U are set to zero to rep- 
resent a complete reflection. The phases (3, 7, 7 consist of 
the dynamical phase and the geometric phase (nonadia- 
batic counterpart to Berry's phase). In physical systems 
(such as the Hubbard model) these phases at the level 
anti-crossings are generally nonzero (and depend on the 
driving force F). 

Expanding the generating function © into a power 
series in z yields the time evolution of the wave func- 
tion (FigEJ) . We can see that \& (distribution on energy 
axis) breaks up, after a transient period, into two parts, 
*(n, t) ~ *cd gc (n-, r)+*bulk(«) t), where VPedgc is a com- 
ponent localized around the boundary (i.e., the ground 
state) while "fbuik is a component traveling into the bulk 
(excited states) with a nearly constant velocity. Interest- 
ingly, the edge state only appears when the phase differ- 
ence, 

= 7-7, 

between the bulk and boundary transfer matrices is 
nonzero (as is the case with physical systems as men- 
tioned). 

What is the nature of the edge state? The two generat- 
ing functions ^f^Q^ n \(z) (cqs.(|3J)) have a common, first- 
order pole in z 2 at z. 2 



polo ~ l-e'Vl^f ' We Can thcn 
obtain, with the Darboux's method[jj|, the asymptotic 
wave function, 4 , e d ge (?T-, t) — ► exp 1 (— ierAt/2h) $ c deeM, 
which obeys Floquet's theorem[ia| with the Floquet 
mode $ and the Floquet quasi-energy per length e/L = 
{F/2tt) argz 2 ole , a function of the electric field F since 
p = cxp(— ttF/F) involves F. We note that e can be 
expressed as e = — ihlim t -too j o\U (t)\^ q) , i.e., the 
asymptotic effective Hamiltonian of the system. An im- 
portant observation here is that the elements of the Flo- 
quet state $ dgo("-) = *(0 L ("), 4> R {n)) with <j> L ' R {n) = 
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edge state propagating slate 



FIG. 2: Time evolution of \ipL,R.{n, r) 2 with the initial con- 
dition i)jl{0, 0) = 1. Here we set p — 0.2 < p c = sin 2 d, and 
only the amplitudes at even n are displayed. 



lim z 2_ >z 2 ^ ^ L ' R (0 — > n;z)(l — z 2 /z 2 olc ), form a geomet- 
ric series, 

|<Mn)| 2 = r n |0 L (O)| 2 , 
\Mn>2)\ 2 = r™" 1 |0 L (O)| 2 , n = 0,2,..., 

r = p/(2-p-2cos6» v /l-p), (7) 

with I (0)| 2 = (1 - r) 2 . Thus $(n) is an edge state ex- 
ponentially localized (on energy axis) around the bound- 
ary n = 0, whose weight is |$ c d go | 2 = X)„(|<Ac(«)| 2 + 
\4>R( n )\ 2 ) = I — r. The size (on energy axis) of the edge 
state is £ = 1/ 1 In r | , which behaves as £ oc F/F in the 
small F regime and diverges like £ oc \p — p c (6)\~ l oc 
F\F - F^e*)!" 1 in_the vicinity of the threshold, p c {6) = 
sin 2 6,F C {0) = -7rF/(21nsin0). When p exceeds p c (9), 
the edge state collapses, and only the component propa- 
gating into the bulk remains. 

Translation to electron systems — Having presented 
the results for the quantum walk, we are now in posi- 
tion to translate them back to the electron system. The 
electric field F causes a production of electron-hole pairs 
through the Landau-Zener tunneling. However, these ex- 
cited charges cannot be accelerated indefinitely, since the 
electron-electron interaction or disorder scatter the mo- 
menta due to the "back scattering" at the level anti- 
crossings. This leads to bifurcations of the amplitude 
in energy space (whose example is depicted in Fig[T^a)). 
When different paths meet, quantum interference induces 
localization of the wave function, and the edge state 
\4>l,r(ji)\ 2 , occurring for 0^0, in the quantum walk 
problem represents an effect of such an interference. So 
the present result is an analytic version of the edge states 
observed numerically in the long-time limit |3J,|J, |7j. 

The momentum and energy expectation values 
(summed over L and R states) for the edge state are 

Jcdgc = £ [jn\Mn)\ 2 +3n\Mn)\ 2 } /|*odgo| 2 , 
n 

£odgc = £[^|^(n)| 2 + ^|^(n)| 2 ]/|$ edge | 2 .(8) 

n 

Here we normalize them by the total amplitude of the 
edge state | ( I?edge('ra)| 2 , which decreases as we increase 
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FIG. 3: Total momentum J e d ge and energy -E c d gc contributed 
by the edge component (eq.JSJ) against the electric field F 
for 9 = 7r/4. The upper figure depicts the weight of the edge 
state. 

the electric field and becomes zero at F = F c (8) (top 
panel of FigEJ. Following [3, 0, 0, El , we choose j£ ,L = 
±jon, Eff> = E^n 2 , n = 0,2,... as the momentum 
and energy of the states, respectively, where jo, Eo are 
units of momentum and energy. We note that the mo- 
mentum measures the chiral asymmetry (between L and 
R) of the distribution. By plugging in eq.lJTJ we ob- 

, • T ■ p(2 — p— 2 cos OyJl — p) J • -1 

tain J c d<?c = 7o n A i /■, a -!v> , an d a similar cx- 
pression for -E G d g c, valid for F < F c {6) (FigE). In 
the weak-field regime, the J dgo is suppressed until the 
Landau-Zener tunneling is activated for F > F when 
Jodgo oc exp(— irF/F) starts to rise. When F is further 
increased to reach F = F c (0), the electric-field induced 
breakdown occurs, which corresponds, in the present pic- 
ture, to the edge-to-propagating transition in the quan- 



tum walk. At the breakdown point the energy expecta- 
tion value diverges as E edge ~ \F — F c (6>)|~ 2 . The J e d ge > 
however, does not diverge but shows a smooth increase. 

Discussions The idea of mapping the non- 
equilibrium problem to a quantum walk may apply to 
wider range of systems having many energy gaps. Here 
we have concentrated on the "edge state" when there is 
randomness (i.e., 6 ^ 0) only at the edge. Randomness 
in the transfer matrices in the excited state, ignored here, 
is expected to enhance the localization to bring |<I> ec z se | 2 
closer to unity. Blatter et.afl have in fact obtained the 
expectation value of the momentum in such a situation, 
where the result resembles J e dge obtained here. However, 
Blatter et.al did not obtain, in the range of electric field 
F they have studied, a dclocalization transition. In fact, 
Cohen et.al, in a random matrix model, found a disap- 
pearance of localized states in strong external fields p. 
If such transitions also occur in electron systems, the 
quantum- walk picture may be used in understanding the 
transition. Indeed, some of the quantum-walk models 
now under study exhibit localization (e.g., |20|). and the 
present mapping may shed lights on the non-equilibrium 
properties of driven quantum systems. 
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